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In perturbative calculations of quantum mechanical path 
integrals in curvilinear coordinates, one encounters Feynman 
diagrams involving multiple temporal integrals over products 
of distributions, which are mathematically undefined. In ad- 
dition, there are terms proportional to powers 5(0) from the 
measure of path integration. We derive simple rules for their 
evaluation from the natural requirement of coordinate inde- 
pendence of the path integrals. 



I. INTRODUCTION 

While quantum mechanical path integrals in curvilin- 
ear coordinates can be defined uniquely and indepen- 
dently of the choice of coordinates within the time-sliced 
formalism [Q, a perturbative definition on a continuous 
time axis poses problems. To exhibit the difficulties, con- 
sider the associated partition function calculated for pe- 
riodic paths on the imaginary-time axis r: 



-A[q] 



where A[q\ is the euclidean action with the general form 
1 



A[q] 



dr 



(2) 



The dots denote r-derivatives, g^i,{q) is a metric, and 
g — det g its determinant. The path integral may for- 
mally be defined perturbatively as follows: The metric 
gfj.:^{q) is expanded around some point Qq in powers of 



. The same thing is done with the poten- 



tial V{q). After this, the action A[q] is separated into a 
free part Ao[qo;Sq] = \g f,i,(.qo)W q" + ^i^^Sqf'Sq" , and 
an interacting part Ant [90; ^q] = A[q] — Ao[qo;Sq]. 

A first problem is encountered in the square root in the 
functional integration measure in (^. Taking it into the 
exponent and expanding it in powers of 6q we define an 
effective action^^ = —^S{0) log[g{qo+Sq)/g{qo)], which 
contains the (5-function at the origin (5(0). It represents 
formally the inverse infinitesimal lattice spacing on the 
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time axis, and is equal to the infinite number 5(0) = 

Jdp/i27T). 

The second problem arises in the expansion of Z in 
powers of the interaction. Performing all Wick contrac- 
tions, Z is expressed as a sum of loop diagrams. There are 
interaction terms involving g^g" which lead to Feynman 
integrals over products of distributions. The diagrams 
contain three types of lines representing the correlation 
functions 



A(r-r')^(g(r)g(r')) = 
drA{r - r') EE (g(T)g(r')) - 
drdr>A{T - r') ^ (g(r)g(r')> - 



(3) 
(4) 
(5) 



The right-hand sides define the line symbols to be used 
in Feynman diagrams for the interaction terms. 
Explictly, the first correlation function reads 



A(r,r') = ^e--l^-^'l. 



(6) 



(j^) The second correlation function (|^) has a discontinuity 



9.A(T,T') = --e(r-r')e-l---'l, 



where 



e(T-r') = -1 



dT"<5(r" - t') 



(7) 



(8) 



is a distribution which vanishes at the origin and is equal 
to ±1 for positive and negative arguments, respectively. 
The third correlation function (|^) contains a (5-function: 



drdr'A{T,T') = S{T-T') 



(9) 



The temporal integrals over products of such distribu- 
tions are undefined 

In this note we define them uniquely by setting up rules 
between these and integrals over products of nonsingular 
correlation functions A(t — r'), plus integrals over pure 
products of (5-functions. These will be defined uniquely 
by the requirement of coordinate invariance of the path 
integral (|l|). 

The internal consistency of these definitions is ensured 
by previous work of the present authors. In Ref. 
we have shown that Feynman integrals in momentum 
space can be uniquely defined as e — > -limits of 1 — e- 
dimensional integrals via an analytic continuation a la 



1 



't Hooft and M. Veltman Q. This definition makes path 
integrals coordinate independent. In Ref. we have 
given rules for calculating the same results directly from 
the Feynman integrals in the 1 — e -dimensional time 
space. 

The present approach has the important advantage 
making superfluous the somewhat tedious analytic con- 
tinuation to 1 — e dimensions. In fact, it does not re- 
quire specifying any regularization scheme. In addition, 
it gives a foundation of a new and general mathematics of 
extending the theory of distributions from a linear space 
to products. 



II. MODEL SYSTEM 

The announced derivation of the identities will be 
based on the requirement of coordinate independence of 
the exactly solvable path integral of a point particle of 
unit mass in a harmonic potential uP'x^ jl^ over a large 
imaginary-time interval /3, 



= Vx{t) e 



The action is 



dT [i2(T) -f wV(t)] 



(10) 



(11) 



A coordinate transformation turns (|T^) into a path in- 
tegral of the type (|^) with a singular perturbation ex- 
pansion. From our work in Rcfs. |^,^ we know that all 
terms in this expansion vanish in dimensional regular- 
ization. Here we shall require the vanishing to find the 
desired identities for integrals over products of distribu- 
tions. 

For simplicity we assume the coordinate transforma- 
tion to preserve the symmetry x <-> —x of the initial 
oscillator, such its power series expansion starts out like 
x{t) = f{q{T)) = q - gq^/3 + g'^aq^/5 - • • • , where g 
is a smallness parameter, and a an extra parameter. We 
shall see that the identities are independent of a, such 
that a will merely serve to check the calculations. The 
transformation changes the partition function (HQ) into 



Z = / Vq{T) e 



-Aj[q]-A[q] 



(12) 



where is A[q] is the transformed action, whereas -A,/ [9] 
an effective action coming from the Jacobian of the co- 
ordinate transformation: 



Aj[q] = -6(0) J dT log 



Sq{T) ■ 



(13) 



The transformed action is decomposed into a free part 



AM = l I dr[<fir)+u;'q'ir)], (14) 
and an interacting part, which reads to second order in 



Aint [q] 



dT { -g 



2e{r)q^r) + ^q\r) 
(l + 2a)^^(r)g4(r)+c.2Q + |)gV) 



(15) 



To the same order in g, the Jacobian action ( p^ ) is 



Aj[q] = -6(0) / dr 



^gq\r) 



(16) 



For (7 = 0, the transformed partition function (|12|) coin- 
cides with (|ol). When expanding Z of Eq. (|l|) m powers 
of g, we obtain sums of Feynman diagrams contributing 
to each order g", which must vanish to ensure coordinate 
invariance. By considering only connected Feynman di- 
agrams, we are dealing directly with the ground state 
energy. 



III. FREE ENERGY DENSITY 

The graphical expansion for the ground state energy 
will be carried here only up to three loops. At any order 
(7", there exist different types Feynman diagrams with 
L = n -|- l,n, and n — 1 number of loops coming from 
the interaction terms (|5|) and (|l^), respectively. The 
diagrams are composed of the three types of lines in (||)- 
(^), and new interaction vertices for each power of g. The 
diagrams coming from the Jacobian action ( |l6|) are easily 
recognized by an accompanying power of (5(0). 

At first order in g, there exists only three diagrams, 
two originated from the interaction (p^), one from the 
Jacobian action ([l^): 



-9 a J 



guj 



00 +9m O- (17) 



At order g'^, we distinguish several contributions. First 
there are two three-loop local diagrams coming from the 
interaction (H), and one two-loop local diagram from the 
Jacobian action (p^): 



15w^ 



1 
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^(0) 00 



(18) 



(19) 



We call a diagram local if it involves no temporal time 
integral. The Jacobian action (16) contributes further 
the nonlocal diagrams: 



g_ 

2! 



25^(0) o - 4(5(0) [ o:> + CO + 2 00 ; 



(20) 
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The remaining diagrams come from the interaction ( |l5| ) 
only. They are either of the three-bubble type, or of the 
watermelon type, each with all possible combinations of 
the three line types The sum of all three-bubbles 

diagrams is 

-|^[ 4 a;x) + 2 oco + 2 cxxj 
+ 8a.^coo+8c.'o:x) +8^' 000 ] , (21) 

while the watermelon-like diagrams contribute 



2! 



(22) 



Since the equal-time expectation value (^(t) ^(t)) van- 
ishes according to Eq. , there are a number of trivially 
vanishing diagrams, which have been omitted. 

In our previous papers all integrals were calcu- 

lated individually in _D = 1 — e dimensions, taking the 
limit e at the end. The results for the integrals en- 
sured that the sum of all Feynman diagrams contributing 
to each order g" vanishes. 



IV. RULES FOR INTEGRALS OVER 
DISTRIBUTIONS 

As a first step in calculating the Feynman integrals 
we express singular time derivatives A(t), A(t) in terms 
of regular correlation functions A(t), plus integrals over 
powers of (5-functions. The tools for this will be partial 
integrations and the inhomogeneous field equation satis- 
fied by the correlation function 



A(r) 



dk 



8{t) + u? 1\(t) , (23) 



Most simply, we have for integrals over products of two 
correlation functions the relation 



dT k^{T)+U?l\^{T) =A(0) 



(24) 



To prove this, we integrate the first term partially, 

drk^ij) = - [ dr A(t)A(t), (25) 



with no boundary term due to the exponential vanishing 
at infinity of all functions involved. Using now the field 
equation (E3h and the property of 5-function that 



dr/(r)<5(T) =/(0), 



(26) 



for any smooth test function /(r), we obtain (24). 

We now turn to singular integrals involving A^ (r) . Us- 
ing the same tools, we obtain, in the same way, the rela- 
tion 



A''{t)+2lu^A''{t)+uj^A''{t) = drS^T). (27) 



dT 



The last integral is undefined. Before fixing its value 
in the next section, we shall derive relations for integrals 
over singular products of four correlation functions. First 
for A(r)A^(T). Using again the field equation (|23|), we 
find 



dT A(t)A3(t) = A3(0) -Lu'^ J dT A\t). (28) 
By a partial integration, the left-hand side becomes 

J dr A(r)A3(r) = "3 ^ dTA\T)A^T), (29) 
leading to 

J drA2(r)A2(T) ^ lA3(0)-iw2y" drA'*(T). (30) 



Invoking once more the field equation (^3|) , we obtain the 
integral 



dTA(r)A'(T)A(T) ^uj' J dTA\T)A\T), (31) 

where we have used A(0) = 0. Due to Eq. (|^), this takes 
the form 

j dr A(r)A2(r)A(r) = icj2A'^(0) - ^uj^ j dT A'^(t). 

(32) 

A further partial integration reduces the integral 



dTA'^iT) 



dT A{t)A^{t)A{t) 



to (^), such that we arrive at the relation 

dTA^{T) = -Uj'^A^{{))+U^ f dTA^{T) 



(33) 



(34) 



We now consider an integral over A2(t)A^ Applying 
again the field equation (ESh , we find the relation 



dTA'^{T)A^{T)= j dTA^T)S^{T) 

- 2lu^A^0) +uj^ f dT A4(t), (35) 



The relations Eqs. (^7|) and ( |35| ) have reduced all in- 
tegrals over singular products of correlation functions to 
regular integrals plus two undefined integrals containing 

(r) . We are now going to show, that the reparametriza- 
tion invariance of path integrals requires the following 
rules for integrals over products of two ^-functions in 
Eqs. (p^) and iM): 



(36) 
(37) 



dT(52(r) =5(0), 
and further 

J dTf{T)6^{T)^f {0)5(0), 

for any smooth test function /(r). 
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V. IMPOSING REPARAMETRIZATION 
INVARIANCE 

To first order in g, the sum of Feynman diagrams ( p^ ) 
must vanisfi: 

O:.) + OO - -5(0) 0=0. (38) 
Tfie analytic form of tliis relation is 

-A(0) + t^2A(o) - 5(0)1 A(0) = 0, (39) 

and the vanishing is a direct consequence of the field 
equation (|3|) for the correlation function at origin. 

At order , the same equation reduces the sum of 
all local diagrams in (^9|) to a finite result plus a term 
proportional to (5(0): 



3Q + a)A(0) + 15(l + ^) 
A2(0) 



-3(a-- )<5(0) 



t^2A(0) 
2 



(40) 



3<5(0)-^c.2A(0) 



Representing right-hand side diagrammatically, we ob- 
tain the identity 



S(19) = 3<S(0) 00-3^' 



(41) 



where E ( [l9| ) denotes the sum of all diagrams in Eq. (19). 
Using the identity (|2j) together with the field equation 
(p3|), we reduce the sum ( pO| ) of all one and two- loop 
bubbles diagrams to terms involving 5(0) and (5^(0): 

-45(0) j dr [a(0) A^ (r) - A(0) A^ (r) + 2ij^A{0)A^ (r 
= 2 A^O) 6(0) + 6^(0) f dT A'' (t). (42) 



(43) 



Hence we find the diagrammatic identity 
1 



2! 



E(20) = 2<5(0) 00 + S'i0)O 



Now, the terms accompanying (5^(0) turn out to be can- 
celed by similar terms coming from the sum of all three- 
loop bubbles diagrams in (^). In fact, the identities ( |2^ ) 
and (p^) lead to 



-4A(0)A(0)A2(t) + 2A2(0)A2(r) 
2A2(0)A2(t) + 8w2a2(o)A2(t) 



8cj2a(0)A(0)A2(t) + 8a;4A2(0)A2(T) 



(44) 



A2(0) -5^(0) / dTA'^ir) 



Thus, we find the diagrammatic identity for all bubbles 
diagrams 

-^^^m~^nM) = - JdrS'ir) CO . (45) 



Finally, the relations (gO|), (|35|), (jS^) and (|3J) reduce 
the sum ( p2[ ) of all watermelon-like diagrams to a finite 
contribution plus the integral involving S'^{t): 



2! 



dr 



A2(T)A2(r) + 4A(T)A2(T)A(r) 

2 

I /I, ,2 A 2/_\ A 2 / ^ 



+ A^(r) + 4a;"A"(T)A"(T) -f ^uj^A^{t) 

O 

= - ^1 dTA^{T)S^{T) + ^CJ^A^{0). 



(46) 



Combining these with all local diagrams (^), we eas- 
ily verify that all finite contributions cancel each other 
leading to the diagrammatic identity 



A2(0). E(19)--E(|l 



3(5(0) -2A-2(0) / dTA^{T)S^{T) 



00 



(47) 



If the singular terms in Eqs. ( p5| ) and (^) are to sum 
up to zero, as required by the coordinate invariance of 
perturbatively defined path integrals, we must have the 
integration rules for the square distribution (^) and (|3^), 
which determined completely the right-hand sides of re- 
lations (|^) and (11). 

The procedure can easily be continued to higher-loop 
diagrams to obtain integrals over any desired products 
of singular correlation functions, and over products of 
(5-functions. 

At no place do we have to specify the value of (5(0) and 
the regularization scheme. There is a perfect cancellation 
of all powers of (5(0) arising from the expansion of the 
Jacobian action, and this is the reason why the so-called 
Veltman rule of setting (5(0) = can be used everywhere 
without problems. 



VI. SUMMARY 

In this note we have set up simple rules for relating 
singular to regular Feynman integrals which avoid the 
explicit calculation of dimensionally regularized integrals 
over products of distributions. These rules follow directly 
from the invariance of perturbatively defined path inte- 
gral under coordinate transformations. Our procedure is 
independent of of regularization prescriptions, using only 
the fact that regularized integrals can be integrated by 
parts. The results are, of course, perfectly compatible 
with those derived before in Refs. |^,|| by dimensional 
regularization. 

Just as in the time-sliced definition of path integrals 
in curved spcae in Ref. j^, there is absolutely no need 
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for extra compensating potential terms found necessary 
in the treatments in Refs. 
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